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Dissociation dynamics of a Bose-Einstein condensate of molecules 
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An unstable condensate of diatomic molecules will coherently disassociate into correlated pairs of 
atoms. This dissociation process exhibits very rich quantum dynamics depending on the quantum 
statistics of the constituent atoms. We show that in the case of bosonic atoms Bose-enhancement 
can lead to stimulated dissociation, whereas, in the case of fermions Pauli-blocking of the available 
states and a build-up of coherence between molecules and atom pairs can give rise to incomplete 
dissociation of the molecules and transient association-dissociation oscillations. 



The ability to create quantum degenerate molecules 
composed of fermionic 0, Q bosonic 0, 0] atoms by 
tuning a molecular level into resonance with the atomic 
states via a Feshbach resonance [1, 0, or by photoas- 
sociation |^, llfl| | has opened up an exciting new area of 
physics for exploration. The case of diatomic bosonic 
molecules coupled to bosonic atoms (6 ^ bb) has been 
shown to undergo coherent association-dissociation oscil- 
lations 1^ and there are predictions of Bose-enhanced 
phenomena in this s yste m that may lead to a, so called, 
superchemistry |3 03 . Il5i | . The case of diatomic 
bosonic molecules coupled to fermionic atoms (6 ^ //) 
has also received a lot of attention lately due to t he p ossi- 
bihty of realizing a BEC-BCS crossover [ll[Tl|l||. For 
positive detuning from resonance (where two-body the- 
ory predicts unstable molecules) molecules can be sta- 
bilized by Pauli-blocking of the atomic states and 
coherent population oscillations have also been predicted 
to occur in this case (20^ 21^ 22] . In these systems quan- 
tum statistics obviously play an important role and it is 
becoming clear that atom-molecule coherence generated 
by their coupling is one of the key elements to under- 
standin g th eir equilibrium and non-equilibrium behavior 



In this paper we consider the production of corre- 
lated pairs of atoms by the spontaneous dissociation of 

Mure Bose-Einstein condensate (BEC) of molecules 
WH . This highly non-equilibrium, spontaneous- 
dissociation regime can be reached by first creating a 
stable BEC of molecules far from the resonance, then 
rapidly tuning through to the other side of the resonance, 
i.e.,j^<0^j^>0, where v is the detuning of the molec- 
ular level from the atomic continuum and depends on the 
applied magnetic field in the case of a Feshbach resonance 
or the detuning of the coupling laser field in the case of 
photoassociation/dissociation. Once the molecule level is 
above the atomic continuum the molecules will become 
unstable and begin to dissociate into atomic pairs. For a 
condensate in a zero momentum state the atoms are cre- 
ated in correlated pairs of equal and opposite momentum 
centered at A: = ^2mv jTi. Correlated pair production by 
this method has been discussed previously in the case of 
bosonic atoms 0, ^3 ■ Here we present a unified treat- 
ment of the dissociation dynamics of a molecular conden- 



sate which includes both the boson (b bb) and fcrmion 
{b ff) case. 

Assuming the momentum spread of the molecular con- 
densate is negligible compared to the mean momentum 
of the emitted atoms, the Hamiltonian of the system can 
be approximated by 

H = hvalao + h'^ujk{bl,^bki + bl^bki) 

k 

+hg^ (^alb^kibk]' + aobl^bl^^'j , (1) 

k 

where ao is the bosonic molecular mode and bk\ and 
bki are the annihilation operators of the atoms and ei- 
ther satisfy bosonic or fermionic commutation relations. 
LUk — hk^ /2m is the dispersion relation of the atoms and 
g is the coupling between the closed channel (molecules) 
and the open channel (free atoms) of the coupled chan- 
nels scattering problem 6]. We have assumed that the 
atoms are in different internal states denoted by t and J, 
but the conclusions can be straight-forwardly applied to 
the case of only one bosonic atomic species. The total 
number of atoms, N = 2aJao -I- X^fcl^-fcT + "fei) con- 
served, where n^g — &;^^&fcs is the number operator for 
the atoms. As atoms of opposite spin and momentum 
are created and destroyed as pairs, the number differ- 
ence, rik] — fi-ki, is also conserved. In the boson case 
the correlation between the atoms created via molecular 
dissociation is analogous to that between photons cre- 
ated in a non-degenerate parametric amplifier (see 
and references within). Quite recently, the pair correla- 
tions between fermionic atoms with opposite momenta 
have been observed in the noise spectrum of photodisso- 
ciated cold molecules [s^l, following the theoretical pro- 
posal of Ref . ■ This pair correlation does not require 
the presence of a molecular condensate as it is simply a 
consequence of the form of the Hamiltonian and arises 
even in an incoherent dissociation process (such as from 
thermal molecules). On the other hand, only coherent 
molecular dissociation, such as from a condensate, can 
give rise to atom-molecule coherence and the related phe- 
nomena of coherent association-dissociation oscilations 
This coherence — characterized by a non-zero value 
of (aJ&_fc|&fcT) — plays an important role in the present 
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work. 

Given that the atoms are created in pairs, we can take 
advantage of a formal mapping between pairs of fermion 
operators and spin- 1/2 Pauli matrices to write the Hamil- 
tonian in a more natural form. This mapping has been 
exploited to determine the phase diagram of the BCS- 
BEC crossover and to predict non-equilibrium atom- 
molecule oscillations ^20, 21, 22]. A similar mapping can 
be made for boson pairs and we treat the two cases in 
parallel. We define new operators by: o'fc_ = &-fc|6fcf, 

crk+ = ^fcT^-fci '^"^ '^'^^ ^ ^("'fcT + ^i-fci T !)• It is 
easy to check that these operators satisfy the commuta- 
tion relations: [crfc±,(Tfc'±] = 0, [a-fcz,crfe'+] = SkM'<^k+, 
[cfczjCTfe'-] = -^fc^fe'Cfc- and 



(2) 



where the upper (lower) sign in Eq. ^ corresponds to 
fermions (bosons). The seemingly insignificant sign dif- 
ference in the commutation relations in the two cases 
leads to completely different dynamics. 

Writing the Hamiltonian Q in terms of these new op- 
erators we have (minus a constant): 



(3) 



and N = 2aQao -l-X]fc(2o'fcz ± 1). In the fermion case, this 
Hamiltonian describes an ensemble of independent two- 
level systems interacting with a single bosonic mode. The 
case of identical two-level systems: ujk — is called the 
Dicke model and is an exactly solvable model that has 
been extensively studied in the quantum optics literature 
(see Ref. [s^l and references within). 

In the special case where the population in each mode 
is small, {{nk] -\- riki) <C 1), throughout the dissocia- 
tion process, due to a large number of available states, 
we can make the approximation that [ak~, Cfe'+l ~ ^k,k' , 
independent of whether the operators describe bosons or 
fermions. In other words, the underlying quantum statis- 
tics of the atoms are unimportant and the Hamiltonian 
describes the coupling of a single mode to a continuum 
of bosonic modes with a quadratic dispersion relation. 
This model has been studied previously in the context 
of an atom laser [s^l- Under the Born-Markov approxi- 
mation, which holds for weak coupling and large v |34| . 
the molecules will experience a rather trivial exponential 
decay and the final frequency distribution of the atom 
pairs will be the standard Lorentzian. 

To proceed further in the general case we note that 
we are interested in the spontaneous dissociation of an 
initially large molecular condensate, so we make a mean- 
field approximation for the molecular mode by replacing 
the operator ao(t) with a c- number a{t) (taken to be 
real without loss of generality) . This approximation will 
break down when the population of molecules approaches 
zero and quantum fluctuations of the molecular mode 
become important. Under the mean-field approximation, 
the atom-molecule coherence reduces to {alb^kibk^) — 




FIG. 1: Example trajectories of Sk{t) for fermions (a) and 
bosons (b). 



a{b-kibk]) = ct{ak~) and the equations of motion for 
the averages (fifes) can be written as: 



dSk 





= ( 2{ujk - v) 




'2{ujk - i 


2ga{t) 



T2ff«(t) 
/ 



Sk, (4) 



where Sk — [{ckx), {cky), (ckz)]* is the column vector of 
averages and we have defined <7kx = (cfc-i- + tTfc-)/2 and 
o'ky = {o'k+ — crfe_)/2i. In addition, a{t) is coupled to 
these variables via number conservation and is given by 



a(t)2^iV/2-^[(afe,(t))±i] 



(5) 



Assuming an initial vacuum state for the atoms, |vac), we 
have (afe,(0)) = (afcy(O)) = and (afc^(O)) = For 
fermions (upper sign) Eqs. Q are the Bloch equations 
describing the dynamics of a two-level system driven 
by the classical field a{t). For each k the motion is 
confined to the surface of the Bloch sphere defined by 
{(Jkx)^ + (cfcy)^ + {(Jkz)'^ = 1/4 and is an expression of 
the underlying Fermi statistics (see Fig.^. On the other 
hand, for bosons (lower sign) the motion is confined to 
the surface of a one-sided three dimensional (3D) hy- 
perbola defined by {crkz)"^ — {^kx)"^ — i'^ky)'^ = 1/4 and 
("■fez) > 1/2, and the population for each k is unbounded 
(see Fig. Pi . 

It is instructive to consider the case when a{t) = ag is 
a constant. In this case, Eq. can be easily solved for 
the above initial state to yield 



1 



(akzit)) ± - = ± 



2f^L 



[l-Ck±{t)], 



(6) 



where r2|_j_ — (.gQ;o)^±(ti'fc— i')^ andCfc+(i) — cos(2r2fe+t) 
and Ck-{t) — cosh(2f2fc_i). Again, the upper (lower) 
signs correspond to fermions (bosons) . In the boson case 
this solution is unstable for 7^ as it leads to exponen- 
tial growth at the rate 2^k- for all > 0. However, 
in the fermion case, the atom population is oscillatory 
about the Lorentzian-shaped mean value: ((7ao)^/2f2|^ 
and, in fact, we can find a self-consistent solution for 
ao by substituting Eq. © back into Eq. (O and assum- 
ing the oscillations eventually dephase for different ujk 
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(this procedure is compared to the numerical solution in 
Fig. |3J). This solution is only valid when there is very 
little molecular decay, but it does indicate that a consis- 
tent solution can be found where the molecule population 
does not completely decay away, and the numerical re- 
sults presented below confirm this. 

For illustrative purposes, we now confine our analysis 
to the case where v is large so that the density of states 
is approximately flat across the region into which the 
molecules tend to decay. In this case, any deviation from 
exponential decay can be directly attributed to the quan- 
tum statistics of the atoms rather than any structure in 
the density of states (c.f. [s^). Taking the continuum 
limit of Eq. lO and evaluating the density of states at v 
we can write 



a{tf^N/2~p[y) d5[{a,{5,t)) ± 



(7) 



where p{v) = V yjir[2{mll\f'l'^ j -k^ is the density of 
states at v for a uniform 3D box of volume V . Since the 
equation of motion for the cr's only depends on (5 = uji^ — v 
we have parameterized them by b instead of fc. As dis- 
cussed above, in the regime where there are ample states 
available to the atoms, the molecules decay exponen- 
tially and the atoms populate the frequencies (5 with a 
Lorentzian distribution of width ^ p{v). A measure of 
the number of available states is therefore given by this 
width multiplied by the density of states: p{v). Moti- 
vated by this, we introduce the dimensionless quantity 
F = N /[gp{v)Y ^ such that the exponential-decay regime 
is given by F <i; 1. It follows that when this condition 
does not hold, we expect the behavior of the dissociation 
process to be altered by the quantum statistics of the 
atoms. 

In Figure |2 and we have plotted the results of a 
numerical solution of Eqs. I@J and Q for F = 4 and 
F = 100, respectively. Already for F = 4 these plots show 
a marked deviation from the usual exponential decay. 
In particular, the molecular population in the fermion 
case does not decay to zero. For F = 100 we see the 
accelerated decay that occurs due to bosonic stimulation 
in the boson case, whereas for the fermion case, a large 
population remains in the molecular state. 

The behavior of the fermionic atoms can be qualita- 
tively understood as follows: The molecular population 
initially undergoes a rapid decay into pairs of fermions. 
However, if F > 1, the states close to resonance, (5 ~ 0, 
become filled and begin to undergo coherent association- 
dissociation oscillations, effectively halting the molecular 
decay. After a few oscillations the molecular population 
settles into a quasi-stationary state, (which we have 
only been able to determine numerically), leaving the 
spin vectors S{5) precessing about the effective "field" 
B{5) = [2gao,Q,25]\ i.e. S{5) = B{5) x S{5). This re- 
sults in the fringes in the distribution over 5 (see insets 
in Figs. 121 and I3J) which become more dense with time. 
The oscillation in the atomic population associated with 
this precession reacts back on the molecular field leading 
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FIG. 2: Molecular population as a function of time for F = 4 
in the case where the decay products are pairs of bosons 
(dashed line) and pairs of fermions (solid line). For com- 
parison the usual exponential decay at the rate 7 — Tvp{u)g^ 
is shown (dotted line). Inset: normalized population distri- 
bution of atoms over 5 = Uk — v at the final time tg — 21.2 
for bosons (dashed line) and fermions (solid line). The dot- 
ted line is a Lorentzian of width 7/4. Only the 5 > case is 
shown as the distribution is symmetric about 5 = 



to an amplitude modulation of the molecular population 
which damps slowly due to the dcphasing of the spins 
with different S. In the final state the spins are com- 
pletely dephased and no net dissociation or association 
can occur. This behavior is reminiscent of the processes 
of optical nutation, where an intense laser pulse excites 
an inhomogeneous media of two- level atoms [ssl ]. Un- 
like the soliton-like oscillations predicted in Refs. [2(tI |. 
the oscillations described here are a transient phenom- 
ena, but have the advantage that they can be created in 
a straight-forward manner, experimentally. The incom- 
plete dissociation for the fermion case was also found in 
Ref. ^] using a stochastic wavefunction approach. 

Due to the assumption that the molecular condensate 
can be described by a mean field the atom pairs with dif- 
ferent k are uncoupled and the Hilbert space of the sys- 
tem can be written as the tensor product H = ^Hk, 
where TCk is the Hilbert space of the pair (fef, —ki). We 
can write the unitary evolution operator as 



Uk{t) ^ exp |rfc(t) 



(8) 

where rfc(t) and 0fe(t) are real time-dependent functions 
and are completely specified by the expectation values 
(o'fcs(i))- Therefore, solving the Eqs. enables us to 
determine not only the expectation values {uks{t)), but 
also the full quantum state of the atoms. In the case 
of bosons, Eq. (jHJ is analogous to the generator of a 
two-photon squeezed vacuum state |2^. In the case of 
fermions, when acting on the vacuum, this evolution op- 
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FIG. 3: Same as Fig. 13 but witii F = 100. Tiie liorizontal 
dotted line shows the steady-state molecular population de- 
termined by a self-consistent solution using Eqs. Q and Q 
(see text) which is a good fit in this case. 



erator yields the state 



\ipkit)) = cosrfc(t) + 6^^6Lfc|e' 



^'^'=(*)sinrfc(t)j |vac), 



(9) 

which has obvious similarities to the BCS state. 

For short times we expect that the dynamics will be 
dominated by the q = condensate of molecules. How- 
ever, the q 7^ molecular modes, collisions and finite 
temperatures that have been neglected from this treat- 
ment are expected to lead to a slow redistribution of the 
atoms and a decay of the correlations Determin- 
ing the relaxation dynamics to the new equilibrium state 
is beyond the scope of the present work, however, we 
can understand some aspects of the relaxation dynam- 
ics by adding phenomenological phase damping terms to 
the equations of motion: d{akj) / dt\rc\ax — ~{<^kj) /T2, 
where j = x, y. This relaxation process leads to a loss 
of the atom-molecule coherence {{<Tk-) 0) and will, 
in turn, destroy the dynamic equilibrium reached in the 



absence of damping. This intuition is confirmed by nu- 
merical simulations which show that for finite T2 the 
molecular population always decays to zero and clearly 
demonstrates the key role played by the atom-molecule 
coherence, in combination with Pauli blocking, in the in- 
complete molecular dissociation in the fermion case. 

Due to their narrow dissociation linewidth, narrow Fes- 
hbach resonances are promising systems to observe the 
effects described here and have the advantage that a mag- 
netic field can be quickly tunned across the resonance into 
the dissociation regime. In fact, the dissociation scheme 
considered here has recently been used to measure the 
width of a number of Feshbach resonances of ®^Rb (38| . 
For the 912 G resonance with a width of AB = 1.3 mG, 
a detuning of hv = kB x l/J-K yields a F > 3 for typical 
densities (n ^ 10^^ cm~'^). Similarly, for the same pa- 
rameters, the narrow 543 G resonance of ^Li which has 
a width of AB = 0.23 G 0], yields a F > 12, demon- 
strating that the regime where quantum statistics play 
a role are well within the reach of current experiments. 
Here we have assumed a uniform 3D system but note 
that the effects of quantum statistics can be significantly 
enhanced in systems of reduced dimensionality or in 
the presence of trappin g po tentials, due to the reduction 
in the density of states |33 ■ 

In summary, we have studied the effects of quantum 
statistics on the dissociation dynamics of a condensate 
of diatomic molecules formed either by two bosonic or 
fermionic atoms. In the former case, the dissociation 
rate is Bose-enhanced; while for the latter, Pauli-blocking 
in combination with the coherence formed between the 
molecules and atom pairs lead to a dynamic equilibrium 
between the molecule and atom populations. Finally, we 
want to point out that we have used a method borrowed 
from quantum optics, which can serve as a powerful tool 
to treat other problems in the coupled atom-molecule 
system, such as the BEC-BCS crossover. 

We thank Randy Hulet for stimulating discussions 
which lead to this work. HP acknowledges support from 
Rice University and Oak Ridge Associated Universities. 



[1] C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Na- 
ture (London) 424, 47 (2003). 

[2] K. E. Strecker, G. B. Partridge, and R. G. Hulet, Phys. 
Rev. Lett. 91, 135 (2003). 

[3] J. CubizoUes, T. Bourdel, S. J. J. M. F. Kokkelmans, 

G. V. Shiyapnikov, and C. Salomon, Phys. Rev. Lett. 
91, 240401 (2003). 

[4] J. Herbig, T. Kraemer, M. Mark, T. Weber, C. Chin, 

H. -C. Nagerl, and R. Grimm, Science 301, 1510 (2003). 
[5] S. Diirr, T. Volz, A. Marte, and G. Rempe, Phys. Rev. 

Lett. 92, 020406 (2004). 
[6] E. Tiesinga, B. V. Verhaar, and H. T. C. Stoof, Phys. 

Rev. A 47, 4114 (1993). 
[7] S. Inouye, M. R. Mathews, J. Stenger, H. J. Miesner, 

D. M. Stamper-Kurn, and W. Ketterle, Nature (London) 

392, 151 (1998). 



[8] E. Timmermans, P. Tommasini, M. Hussein, and A. Ker- 

men, Phys. Rep. 315, 199 (1999). 
[9] P. O. Fedichev, Y. Kagan, G. V. Shiyapnikov, and 

J. T. M. Wahaven, Phys. Rev. Lett. 77, 2913 (1996). 
[10] F. K. Fatemi, K. M. Jones, and P. D. Lett, Phys. Rev. 

Lett. 85, 4462 (2000). 
[11] E. A. Donley, N. R. Claussen, S. T. Thompson, and C. E. 

Wieman, Nature (London) 417, 529 (2002). 
[12] D. J. Heinzen, R. Wynar, P. D. Drummond, and K. V. 

Kheruntsyan, Phys. Rev. Lett. 84, 5029 (2000). 
[13] M. G. Moore and A. Vardi, Phys. Rev. Lett. 88, 160402 

(2002). 

[14] K. V. Kheruntsyan and P. D. Drummond, Phys. Rev. A 

66, R031602 (2002). 
[15] A. Vardi and M. G. Moore, Phys. Rev. Lett. 89, 090403 

(2002). 



5 



[16] A. J. Leggctt, Modem Trends in the Therory of Con- 
densed Matter ( Springer- Vcrlag, 1980), edited by A. 
Pekalski and R. Przystawa. 

[17] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 
92, 040403 (2004). 

[18] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 
Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. 
Lett. 92, 120403 (2004). 

[19] G. M. Falco and H. T. C. Stoof, Phys. Rev. Lett. 92, 
0130401 (2004). 

[20] A. V. Andreev, V. Gurarie, and L. Radzihovsky, Phys. 
Rev. Lett. 93, 130402 (2004). 

[21] R. A. Barankov and L. S. Levitov, Phys. Rev. Lett. 93, 
130403 (2004). 

[22] C. P. Search, W. Zhang, and P. Meystre, Phys. Rev. A 

91, 190401 (2004). 
[23] M. Holland, S. J. J. M. F. Kokkelmans, M. L. Chiofalo, 

and R. Walser, Phys. Rev. Lett. 87, 120406 (2001). 
[24] M. Mackie, K.-A. Suominen, and J. Javanainen, Phys. 

Re. Lett. 89, 180403 (2002). 
[25] M. Greiner, C. A. Regal, and D. S. Jin, Nature (London) 

426, 537 (2003). 
[26] S. Jochim, M. Bartenstein, A. Altmeyer, G. Hendl, 

C. Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. 

Lett. 91, 240402 (2003). 
[27] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 

Raupach, S. Gupta, Z. Hadzibabic, and W. Ketterle, 



Phys. Rev. Lett. 91, 250401 (2003). 
[28] T. Bourdel, L. Kaykovich, J. Cubizolles, J. Zhang, 
F. Chevy, M. Teiehmann, L. Tarruell, S. J. J. M. F. 
Kokkelmans, and C. Salomon, Phys. Rev. Lett. 93, 
050401 (2004). 

[29] D. F. Walls and G. J. Milburn, Quantum Optics 
(Springer- Verlag, New York, 1994). 

[30] M. Greiner, C. A. Regal, J. T. Stewart, and D. S. Jin, 
Phys. Rev. Lett. 94, 110401 (2005). 

[31] E. Altinan, E. Dernier, and M. D. Lukin, Phys. Rev. A 
70, 013603 (2004). 

[32] S. M. Chumakov and M. Kozierowski, Quantum Semi- 
class. Opt. 8, 775 (1996). 

[33] J. J. Hope, Phys. Rev. A 55, R2531 (1997). 

[34] M. W. Jack, M. J. CoUett, and D. F. Walls, Phys. Rev. 
A 59, 2962 (1999). 

[35] L. Allen and J. H. Eberly, Optical Resonance and Two- 
Level Atoms (John Wiley and Sons, New York, 1975). 

[36] J. F. Corney and P. D. Drummond, Phys. Rev. Lett. 93, 
260401 (2004). 

[37] V. M. Galitski, Phys. Rev. A 72, 013612 (2005). 

[38] S. Diirr, T. Volz, and G. Rempe, Phys. Rev. A 70, 
031601(R) (2004). 

[39] C. P. Search, H. Pu, W. Zhang, B. P. Anderson, and 
P. Meystre, Phys. Rev. A 65, 063616 (2002). 



